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摘要 :【 目 的) 鉴定 中 华 按 蚁 Anopheles sinensis , F] VG 3E. 32: 3X. Anopheles gambiae 35 RAF 33. Aedes aegypti 
7a Sk £8 Culex quinquefasciatus 4 p ERE FRIDA E 89 AF RRA IR8a 与 IR25a 以 及 代表 
性 的 iGluRs 基因 和 IRs 基因 ,研究 这 些 基因 的 特征 及 其 系统 发 育 关系 ,确定 IRSa 与 IR25a 基因 在 
离子 谷 氨 酸 受 体 基因 中 的 分 类 地 位 。【 方 法】 以 黑 腹 果 晶 Drosophila melanogaster 的 iGluRs 和 IRs 
基因 序列 作为 询问 序列 ,用 BLASTP 和 BLASTN Zr ;k 2-3 34 A P AE di, EG E dC IRURE C fe 
致 倦 库 蚁 的 基因 组 数据 库 , 以 鉴定 这 4 个 种 内 的 icluRs fe IRs 基因 序列 ;用 生物 信息 学 方法 比较 分 
析 鉴 定 获得 的 IRSa 与 IR25a 及 代表 性 的 iGluRs 和 IRs 基因 和 氨基酸 序列 的 特征 ;用 最 大 似 然 法 、 贝 
叶 斯 法 和 最 大 简约 法 构建 并 讨论 这 些 基 因 和 氨基 酸 序 列 的 系统 发 育 关 系 ; 用 PAML 软件 计算 这 些 基 
的 dwds(@) 值 并 分 析 其 选择 压力 。 【结果 】 通 过 对 中 华 按 明 、 冈 比 亚 按 虹 、 埃 及 伊 虹 和 致 倦 库 玉 
中 离子 谷 氨 酸 受 体 基因 进行 生物 信息 学 鉴定 ,获得 了 所 有 IRSa 和 IR25a 基因 ,24 个 代表 性 的 
iGluRs 基因 和 30 个 代表 性 的 IRs 基因 的 序列 。 特 征 比较 和 系统 发 育 关系 分 析 结 果 表 明 ,在 4 种 蚁 
忠和 黑 腹 果 蝇 中 IRSa 与 IR25a 基因 的 氨基 酸 序列 长 度 与 iGluRs 的 更 接近 ,其 平均 值 远 大 于 其 他 
IRs 45 ; IR25a 和 IR8a 有 与 传统 iGluRs 相同 的 氨基 酸 末 端 区 域 (amino-terminal domain, ATD) , 虽然 
IR8a 的 序列 不 太保 守 , 而 IRs 则 没有 ATD;IR8a 与 IR25a 序列 有 与 传统 iGluRs 一 样 保守 的 配 体 结 
AR SI 和 S2 ,而 JIRs 没 有 。 所 有 这 些 基 因 的 o 值 远 小 于 1 ,意味 着 这 些 基因 在 进化 过 程 中 都 经 历 
了 纯化 选择 ,而 IR8a 与 IR25a 的 o@ 值 更 接近 iGluRs。 此 外 ,IR8a 与 IR25a 还 有 与 non-NMDA 相似 
的 特异 性 位 点 ,而 这 些 位 点 在 IRs 不 存在 。 再 者 ,IR8a 和 IR25a 与 non-NMDA 聚 在 一 枝 ,形成 一 个 
姐妹 群 , 群 内 序列 间 的 遗传 距离 小 于 本 研究 中 任何 其 他 序列 间 的 距离 。 根 据 这 些 研究 结果 ,我 们 将 
IR8a 5 IR25a 分 类 进入 传统 iGluRs 中 的 non-NMDA 类 型 ,并 命名 为 Putative 受 体 。【 讨 论 ] 本 研究 
13i Y 3x x iGluRs 和 IRs 基因 的 信息 框架 ,确定 了 IRSa 与 IR25a 基因 的 分 类 地 位 ,对 进一步 的 功 
能 研究 具有 重要 意义 。 

关键 词 : Eo. iGluRs; IRs; JR8a; IR25a; 全 基因 组 鉴定 ; 序列 特征 ; 系统 发 育 
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Abstract: [ Aim] To identify the ionotropic receptor genes /R8a and IR25a, and the representative 
iGluRs and IRs genes in four vector mosquito species of medical importance, Anopheles sinensis, 


Anopheles gambiae, Aedes aegypti and Culex quinquefasciatus, and to investigate the characteristics and 
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phylogenetic relationships of these genes, so as to determine the classification position of the /R8a and 
IR25a genes in the group of ionotropic glutamate receptor genes. [Methods] By using BLASTP and 
BLASTN searching against the genomes of these four mosquito species with the iCluRs and IRs gene 
sequences of Drosophila melanogaster as query sequences, the sequences of iGluRs and IRs genes of 
these four mosquito species were identified. The characteristics of the amino acid sequences of IR8a and 
IR25a and the representative iGluRs and IRs genes identified were predicted and comparatively analyzed 
by bioinformatics methods. The phylogenetic relationships of these genes in the five species were deduced 
and discussed by using maximum likelihood, Bayesian and maximum parsimony methods based on amino 
acid sequences, and the selection pressures of these genes were analyzed with the dy⁄/ds (œ) values 
calculated by using PAML software. [Results] The sequences of /R8a and IR25a, 24 representative 
traditional iGluRs genes and 30 representative IRs genes in mosquitoes An. sinensis, An. gambiae, An. 
aegypti and Cx. quinquefasciatus were identified by bioinformatics methods. The comparative analyses of 
characteristics and phylogenetic relationships of these genes showed that the amino acid sequence lengths 
of IR8a and IR25a in the four mosquito species and D. melanogaster are closer to those of iGluRs, and 
their mean values are significantly longer than those of other IRs. IR25a has the amino-terminal domain 
( ATD) that is the characteristic of traditional iGluRs, and IR8a also has a region similar to ATD, 
although it is not conservative, while IRs have no ATD. IR8a and IR25 also have the conserved S1 and 
S2 ligand-binding domain ( LBD) , the characteristics of iGluRs, but IRs do not have it. The d,/ds( œ) 
values of all these genes were far less than 1, suggesting the purifying selection in evolution, while the 
d/d; (œw) values of IR8a and IR25a were closer to those of iGluRs. In addition, like non-NMDA, IR8a 
and IR25a have specific sites which do not exit in IRs. IR8a and IR25a cluster together with non-NMDA 
to form a sister group, with the genetic distance between the sequences within the group less than that of 
any other pairs of sequences in this study. Based on these results, IR8a and IR25a were classified into the 
non-NMDA receptors in iGluRs, and named as Putative receptors. [Conclusion] This study developed the 
information frame of iGluRs and IRs genes of mosquito species and determined the classification position of 
IRS8a and IR25a genes, which is of significance for further functional study of these genes. 

Key iGluRs; IRs; IR8a; IR25a; 


characteristics; phylogeny 


words: Culicidae; genome-wide identification; sequence 


离子 谷 氮 酸 受 体 (ionotropic glutamate receptors , 2004; Ashby et al., 2008) 。 











iGluRs) 是 真 核 生物 和 原核 生物 中 一 类 保守 的 配 体 
结合 离子 通道 ,在 疹 椎 动物 神经 系统 突 触 间 的 信号 
传递 中 作为 兴奋 的 神经 递 质 谷 氨 酸 的 受 体 ,通过 调 
节 阳 离子 流 跨 质 膜 来 应 答 细胞 外 谷 氨 酸 和 相关 配 体 
的 结合 ,细胞 可 以 利用 其 来 调节 内 部 生理 以 应 对 外 
部 环境 信号 (Mayer and Armstrong, 2004) , iGluRs 
在 药理 学 和 分 子 特 征 上 被 分 为 3 个 亚 类 ,根据 它们 
的 主要 激动 剂 将 其 分 别 命名 为 :a- 氮 基 -3 羟基 -5 甲 
IEA 异 恶 唑 受 体 (AMPA receptors ) , 4T. $E 24 2 52 1% 
(Kainate receptors ) 和 N-FH JE-D-X 4€ A We Az pk 























离子 受 体 ( ionotropic receptors, IRs) Z& J Jie 
iGluRs 中 新 发 现 的 一 类 基因 (Benton et al., 2009) , 
与 嗅觉 相关 (Vosshall and Stocker, 2007; Masse et 
al., 2009) ,在 既 不 表达 嗅觉 受 体 (olfactory receptor, 
OR) 基因 也 不 表达 味觉 受 体 ( gustatory receptor, 
GR) 基因 的 触角 腔 锥 形 神经 元 中 表达 ( Vosshall and 
Stocker, 2007 ; Montell, 2009) ,因此 被 认为 是 一 类 
新 的 化 学 感应 受 体 基因 ,在 检测 外 界 环境 中 的 挥发 
性 化 学 物质 以 及 细胞 间 化 学 信号 的 过 程 中 有 广泛 的 
功能 (Benton，2007 ) 。 传 统 的 iGluRs 和 IRs 已 在 许 











(NMDA receptors) 。AMPA 受 体 主要 在 大 部 分 中 枢 
神经 系统 突 触 内 调节 快速 兴奋 性 突 触 传导 , 而 
Kainate 受 体 在 这 个 过 程 中 只 有 微弱 的 调节 作用 , 它 
们 通常 被 合 称 为 non-NMDA 受 体 。NMDA 受 体 则 在 
突 触 的 可 塑性 (synaptic plasticity ) 中 发 挥 功能 , 需 

谷 氨 酸 和 甘氨酸 两 种 激活 剂 来 激活 ( Heidelberg, 





多 动物 物种 中 被 鉴定 出 来 ( Tikhonov and Magazanik , 
2009) 。 在 黑 腹 果 晶 Drosophila melanogaster、 风 比 亚 
按 蚊 Anopheles gambiae RNAF IX Aedes aegypti HS 
f FEX. Culex quinquefasciatus 基因 组 中 ,分 别 有 12, 
9,8 和 6 个 iGluRs 基因 及 66, 46, 94 和 69 个 IRs 
基因 (Liu et al., 2010; Croset et al., 2010) 。 在 所 有 
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已 鉴定 IRs 基因 的 昆虫 中 ,都 有 IR8a 与 IR25a 基因 存 
在 ,它们 在 触角 的 大 部 分 位 置 单独 或 与 别 的 IR 基因 
一 起 表达 ( Benton et al., 2009) ,被 归 类 为 IRs 基因 。 
但 是 由 于 IR8a 与 IR25a 基因 在 序列 特征 .保守 性 和 
系统 发 育 关系 等 方面 的 特异 性 ,它们 在 分 类 家 族 的 归 
属 上 一 直 是 一 个 未 解决 的 问题 (Croset et al., 2010) 。 

首先 在 氮 基 酸 序列 上 ,传统 iGluRs ( NMDA, 
non-NMDA ) 和 IRs 都 包含 一 个 胞 外 N- 末 端 ,一 个 被 
离子 通道 孔 P 隔 开 的 由 Sl 和 S2 构成 的 配 体 结合 区 
(ligand binding domain, LBD ) 
Armstrong, 2004; Jin et al., 2009) ,3 个 跨 膜 区 (M1 — 
M3) 和 一 个 胞 内 的 C- 末 端 ( Mayer et al., 2006)。 其 
中 离子 通道 孔 P 和 跨 膜 区 MI — M3 组 成 离子 通道 
结构 域 Lig Chan ( PF00060) (Kuner et al., 2003) 。 
此 外 ,传统 的 iGluRs TE N 末端 还 有 一 个 氨基 酸 末 端 
区 域 ( amino-terminal domain, ATD ) ( PF01094 ) , ffi 
责 通道 的 组 装 以 及 结合 辅助 因子 。 黑 腹 果 晶 的 
IR8a 和 IR25a 也 有 与 ATD 相似 的 序列 ,只 是 更 分 散 
( Croset et al., 2010) ,然而 其 他 的 IRs (不 包括 IRSa 
和 IR25a, 下 同 ) 则 没有 这 个 序列 。 其 次 , 传统 
iGluRs (不 包括 IRs) 的 配 体 结合 位 点 保守 性 强 , 它 
们 在 Sl 区 都 有 与 谷 氨 酸 盐 的 a 羧基 直接 结合 的 精 
ARR) ,S2 区 都 有 与 谷 氮 酸 盐 的 y 羧基 直接 结合 
£15] 2) 24 82 CT) 和 与 谷 氮 酸 盐 的 a 氨基 酸 直接 结合 的 
天 冬 氮 酸 (D) 或 谷 氮 酸 (上 )(Naur et al., 2007) 。 黑 
JE wA. uw g Acyrthosiphon pisum, JR UG XX 
Tribolium castaneum ,JK BT Y Spodoptera littoralis 和 
ZÆ Bombyx mori 中 的 IR8a 同时 有 3 个 配 体 结合 位 
点 ,IR25a 也 有 其 中 两 个 (T,，D/E)。 其 他 IRs 的 配 
体 结合 区 则 非常 分 散 , 大 部 分 都 没有 或 者 只 有 其 中 
一 个 位 点 (Benton et al., 2009; 郭 梅 , 2014) ,可 能 结 
合 其 他 类 型 的 配 体 行 使 不 同 的 功能 。 再 者 ,在 黑 腹 
果 蝇 中 基于 氨基 酸 序列 构建 的 iCluRs，IRs，IR8a 
和 IR25a 的 系统 发 育 上 ,IR8a 与 IR25a 不 聚集 在 IR 
家 族 支 上 ,而 在 iGluRs 3c E ( Bengtsson et al., 2012; 
Liu et al., 2012) ,进化 速率 也 与 iGluRs 更 接近 , 远 
慢 于 其 他 的 IRs 基因 。 这 种 情况 在 坚 虫 和 棉铃 虫 
Helicoverpa armigera 等 物种 的 系统 发 育 关 系 中 也 有 
发 现 (Liu et al., 2012; Glaser et al., 2013) 。 

本 研究 基于 重庆 师范 大 学 中 华 按 蚊 基因 组 和 转 
录 组 测序 数据 (Chen et al., 2014) 及 从 Vectorbase 下 
FREY b] EE ME Tz .埃及 伊 蚊 和 致 倦 库 蚊 的 基因 组 数 
据 , 以 黑 腹 果 蝇 iGluRs 和 IRs 基因 序列 作为 询问 序 
列 ,采用 BLASTN 和 BLASTP 方法 ,从 信息 学 上 鉴定 
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了 这 4 种 重要 医学 媒介 蚊虫 的 IRSa 与 IR25a 基因 
以 及 代表 性 的 iGluRs 和 IRs 基因 ,并 进行 一 系列 特 
征 分 析 。 通 过 将 IRSa 和 IR25a 与 iGluRs 和 IRs 从 
氨基 酸 序 列 长 度 、 结 构 域 .保守 性 替换 率 、 特 异性 位 
点 和 系统 发 育 关系 6 个 方面 进行 比较 分 析 , 确认 了 
IR8a 与 IR25a 的 特殊 性 ,并 将 它们 归 为 iGluRs 中 的 
non-NMDA 受 体 类 型 ,命名 为 Putative 受 体 。 此 研究 
丰富 了 蚊虫 的 iGluRs 基因 数据 ,为 进一步 研究 IRSa 
与 IR25a 基因 提供 了 信息 基础 。 




















1 材料 与 方法 


1.1 数据 来 源 

中 华 按 蚊 基因 组 数据 包括 基因 组 序列 和 注释 以 
及 预测 的 基因 集 序列 ,由 重庆 师范 大 学 昆虫 和 分 子 
生物 学 研究 所 在 华 大 基因 完成 ,组 装 覆 盖 度 达到 了 
转录 组 测序 产生 的 蛋白 编码 ESTs 的 98. 3596 ( Chen 
et al., 2014)。 冈 比 亚 按 蚁 、 埃 及 伊 改 和 致 倦 库 蚊 的 
基因 组 数据 库 和 黑 腹 果 蝇 的 iGluRs ,IR8a 和 IR25a 
以 及 IRs 基因 序列 下 载 自 Vectorbase (https://www. 
vectorbase. org/ ) o 
1.2 基因 的 生物 信息 学 鉴定 

首先 ,以 黑 腹 果 蝇 iGluRs 和 IRs 基因 作为 查询 
序列 ,对 上 述 4 种 蚊虫 的 基因 组 数据 库 进行 
BLASTP 和 BLASTN 同 源 性 搜索 ,E-value < 1 x 10 7 
作为 闵 值 ,将 得 到 的 候选 基因 在 Softberry 中 预测 其 
转录 组 和 和 蛋白 序列 后 ,于 NCBI 中 再 进行 BLASP 搜 
索 ,查看 其 结构 域 以 及 与 查询 序列 的 相似 度 ,鉴定 得 
到 正确 的 基因 。 
1.3 基因 特征 分 析 

使 用 BLAST TE (http://www. ncbi. nlm. nih. 
gov/ BLAST/) 进行 序列 相似 性 搜索 ;使 用 NCBI Blast 
和 Interproscan( http :// www. ebi. ac. uk/interpro/ ) 对 
氨基 酸 序列 进行 结构 域 预 测 和 亚 细 胞 定位 。 使 用 
Clustal W 软件 (Thompson et al., 2002) 对 4 种 蚊虫 
和 黑 腹 果 蝇 的 代表 性 iGluRs, IR8a, IR25a 和 IRs 的 
氨基 酸 序 列 进行 多 重 比 对 。Gap Opening Penalty 参 
数 设 为 10, Gap Extension Penalty iX Jy 0.2; Gap 
Separation Distance ix 7j 4; Delay Divergent Cutoff 
(96 ) 设 为 30% 。 和 蛋白 质 的 三 级 结构 ,要 首先 通过 
PSI-BLAST 在 PDB ( Protein Data Bank ) 中 搜索 得 到 
氮 基 酸 序 列 一 致 性 最 高 (大 于 80% ) , 且 具 有 相应 结 
构 域 的 基因 ,用 其 作为 同 源 建 模 的 模板 , 再 采用 
SWISS-MODEL ( http ;//swissmodel. expasy. org/) ilt 
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行 同 源 建 模 和 3D 结构 预测 ( 唐 尧 等 , 2014) 。 
1.4 系统 发 育 分 析 

Ti rp edt, px pe ME 2 SC, SUP Poe IR UD 
和 黑 腹 果 蝇 5 种 昆虫 的 IR8a 和 IR25a (3510 4) K 
代表 性 的 iGluRs(24 个 ) 和 代表 性 的 IRs (30 个 ) 24 
基 酸 序列 ,用 Clustal W 对 其 进行 多 重 比 对 ,分 别 用 
最 大 似 然 法 (maximum likelihood, ML) 、 贝 叶 斯 法 
(Bayesian inference, BI) 和 最 大 简约 法 (maximum 
parsimony, MP) 构建 系统 发 育 树 。 其 中 ,用 Mega 
5.0 软件 构建 ML 和 MP 系统 发 育 树 ,用 Find Best 
Protein Model 进行 模型 预测 ,ML 树 用 WAG 模型 。 
使 用 1 000 次 重复 计算 系统 发 育 树 上 的 bootstrap 值 
( 自 展 分 析 值 )。 使 用 PAUP 软件 构建 贝 叶 斯 系统 
发 育 树 ,用 Modeltest3.7 软件 分 析 模 型 ,BI 树 也 用 
WAG 模型 。 大 于 5096 的 自 展 值 标 于 系统 发 育 树 相 
应 分 文 上 ,讨论 其 系统 发 育 关系 。 
1.5 基因 替换 率 分 析 

使 用 PAML( Yang, 1997) 软件 中 的 最 大 似 然 法 
计算 基因 的 dv/ds(。w) 值 。 以 中 华 按 蚊 的 序列 作为 
对 照 , 用 MEGA 软件 中 的 MUSCLE codon 分 别 将 其 
与 冈比亚 按 蚊 , 致 倦 库 蚊 和 埃及 伊 蚊 直系 同 源 基因 
的 编码 区 (CDs) 进行 比 对 ,手动 删除 其 中 的 终止 子 ， 
存 为 FASTA 格式 。 将 FASTA 格式 转换 成 PAML 格 
式 后 ,输入 PAML 软件 ,于 PAML 的 yn00 中 ,用 MO 
模型 计算 非 同 义 突 变频 率 (dy)、 同 义 突 变频 率 
(ds) 和 d/ds (w) RE, o1 为 纯化 选择 ,@ 关 1 为 
正 向 选择 ,w =1 为 中 性 选择 ( Swanson and Vacquier, 











2002; Swanson et al., 2003) 。 本 研究 计算 了 3 种 蚊 
虫 中 的 IRSa 和 IR25a 基因 及 代表 性 iGluRs 和 IRs 
基因 的 o 值 , 比 对 分 析 所 受到 的 选择 压力 和 保 
守 性 。 


2 结果 


2.1 基因 的 生物 信息 学 鉴定 

本 研究 通过 对 中 华 按 蚊 px EGAV P c RKI 
和 致 倦 库 蚊 中 的 离子 谷 氨 酸 受 体 基因 进行 生物 信息 
学 鉴定 ,获得 了 所 有 IR8a 和 IR25a 基因 ,以 及 24 个 
代表 性 的 iGluRs 基因 和 30 个 代表 性 的 IRs 基因 的 
序列 信息 。 其 中 致 倦 库 蚊 CglR21a, CqIR25a 和 
CqIR68a 是 非 全 长 序列 ,因此 没有 包括 在 后 续 3 类 
基因 氨基 酸 序列 长 度 和 dvds(o) 值 的 比较 计算 中 。 
2.2 IR8a 5E IR25a 基因 的 结构 

基于 4 种 蚊虫 和 黑 腹 果 晶 IRSa 和 IR25a 基因 
组 序列 和 转录 组 序列 ,在 Gene Structure Display 
Server 在 线 网 站 中 预测 得 到 如 图 1 所 示 的 基因 结 
构 。IR8a 基因 有 3 ~5 个 外 显 子 ( 图 1: A),IR25a 
基因 有 8 ~ 10 个 外 显 子 (图 1: B) 。 对 于 IR8a 和 
IR25a 两 个 基因 的 内 含 子 ,长 度 在 50 ~ 100 bp 的 个 
数 占 58.18% ,长 度 在 100 ~ 200 bp 的 个 数 占 
21.8296 (图 2);0, 1 和 2 相位 内 含 子 所 占 比 例 分 
别 为 30.91% , 34. 55% FI 34. 5596 , AaIR8a 基因 的 
结构 与 其 他 基因 差别 较 大 ,有 一 个 大 于 1 000 bp 的 
内 含 子 (图 1: A). 
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ED 外 显 子 Exon — NÉT Intron 0, 1, 2: 内 含 子 相位 Intron phase 
图 1 5 种 昆虫 IR8a 与 IR25a 的 基因 结构 
Fig. 1 Gene structure of /R6a and [R25a in five insect species 


As; 中 华 按 蚊 Anopheles sinensis; Ag: 冈比亚 按 蚊 Anopheles gambiae; Aa; 埃及 伊 蚊 Aedes aegypti; Cq: 致 倦 库 蚊 Culex quinquefasciatus; Dm; 黑 腹 
TUB Drosophila melanogaster. 外 显 子 用 蓝 色 的 圆 角 和 矩形 表示 ,内 含 子 用 黑色 的 线段 表示 , 线 上 的 数字 表示 各 内 含 子 的 相位 。Exons are marked 


with blue rounded rectangles, introns are marked with black lines, and the numbers above the black lines indicate intron phase. 
































4 期 王 婷 婷 等 : 四 种 重要 医学 媒介 蚊虫 离子 受 体 基 因 IR8a 和 IR25a 的 特征 及 分 类 地 位 383 






















58.1896 


tA 
© 


内 含 子 频率 (%) 
Frequency of introns 
人 
5 


: 20 
27% 9.0996 


1-50 50-100 100-200 200-500 500-1000 21000 





内 含 子 长 度 范围 (bp) 


Range of intron size 





图 2 5 种 昆虫 IR8a 与 IR25a 基因 
不 同 长 度 内 含 子 的 频率 


Fig. 2 Frequencies of different lengths of introns of /R8a 





and /R25a in five insect species 


2.3 iGluRs 和 IRs 氨基 酸 序列 长 度 比较 

4 种 蚊虫 和 黑 腹 果 蝇 iGluRs 的 氨基 酸 序列 长 度 
在 632 ~1 055 之 间 ( 平 均值 =946 a. a. ) ; IR8a 与 
IR25a 的 氨基 酸 序列 长 度 在 887 ~ 1 116 之 间 ( 平 均 
值 =940 a. a. ) ,而 IRs 的 氨基 酸 序列 较 短 ,在 481 ~ 
1 047 之 间 ( 平 均值 =667 a. a. ) ,表明 IR8a 和 IR25a 
基因 在 氨基 酸 序 列 长 度 上 与 iGluRs 的 平均 值 更 为 
接近 (图 6: A). 
2.4 iGluRs 和 IRs 氨基 酸 序 列 结 构 域 . 配 体 结合 
位 点 和 特异 性 位 点 比较 
2.4.1 结构 域 比较 : 传统 iGluRs 和 IR25a( 不 包括 
非 全 长 CqIR25a) 氨基 酸 序列 都 有 和 氨基酸 末 端 区 域 






































(ATD) . 配 体 结合 区 (SI1，S2) 和 离子 通道 区 Lig_ 


Chan 3 个 结构 域 ; IR8a 有 配 体 结合 区 和 Lig_Chan 
结构 域 ;而 IRs 则 只 有 Lig_Chan 结构 域 (图 3: B)。 

2.4.2. 配 体 结合 位 点 比较 : 传统 的 iGluRs 与 IR8a 
氨基 酸 序列 在 配 体 结合 区 都 有 直接 与 谷 氮 酸 配 体 结 
合 的 3 个 氨基 酸 位 点 : 精 氨 酸 (R)、 苏 氨 酸 (T) 、 天 
ZAR) EAT 2482 (E), IR25a FFARR) 和 
天 冬 氨 酸 (D) 两 个 氨基 酸 位 点 ,IRs 无 这 3 个 位 点 
(图 3: A)。 这 说 明 IR8a 与 IR25a 的 配 体 结合 位 点 
保守 性 强 , 不 同 于 配 体 结合 区 极度 分 散 的 IRs 
(Glaser et al., 2013; Liu et al., 2014) , 

2.4.3 特异 性 位 点 比较 : 传统 iGluRs 的 氨基 酸 序 
列 有 一 些 特异 性 位 点 (special sites) ,在 动物 神经 元 
突 触 内 以 及 谷 氮 酸 激活 钟离 子 通道 的 过 程 中 发 挥 作 
用 。 这 些 特 异性 位 点 有 两 个 主要 的 特征 :多 肽 结合 
位 点 ( peptide binding site) ( Brockie and Maricq, 
2003) FU — XE pk 7 qi (dimer interface) ( Madden, 


















































2002) ,具有 这 些 位 点 的 氨基 酸 序 列 属于 胞 质 结合 
蛋白 2 超 家 族 的 成 员 。Kainate 和 AMPA 受 体 的 特异 
性 位 点 分 别 是 PBP2_iGluRs_Kainate ( cd13714 ) 和 
PBP2_iCluRs_AMPA (cd13715) ,NMDA 受 体 的 特异 
性 位 点 是 PBP2_iGluRs_NMDA ,根据 组 成 亚 基 的 不 同 
又 分 为 NMDA _ Nrl/NI2/NI3. ( cd13719/cd06378/ 
cd13720) 3 种 。 我 们 在 IR8a 与 IR25a 氨基 酸 序列 中 
发 现 了 PBP2 iGluRs putative (cd13717 ) 特异 性 位 点 ， 
而 在 IRs 氨基 酸 序列 中 没有 发 现任 何 特异 性 位 点 。 

通过 氨基 酸 序列 多 重 比 对 发 现 ,上 述 所 有 这 些 
寺 异 性 位 点 都 出 现在 配 体 结合 区 SI 和 S2 中 ,IR8a 
和 IR25a 单独 与 传统 non-NMDA 受 体 相似 的 位 点 有 
12 个 (序列 上 方 红 色 三 角形 标记 的 位 点 ) ,同时 与 传 
统 non-NMDA 和 NMDA 受 体 相似 的 位 点 有 11 个 
(序列 上 方 蓝 色 圆 形 标记 的 位 点 ) ,但 是 没有 单独 与 
NMDA 受 体 相似 的 位 点 (图 3: A)。 另 外 我 们 在 胞 
质 结 合 蛋白 2 超 家 族 的 分 层 CD 树 (http://www. 
ncbi. nlm. nih. gov/Structure/cdd/eddsrv. cgi) 中 可 以 
看 到 , 具有 PBP2 _iGluRs Kainate, PBP2 _iCluRs _ 
AMPA 和 PBP2. iGluRs, putative 特异 性 位 点 的 氨基 
酸 序列 都 归 类 为 non-NMDA 受 体 (Brockie and 
Marieq, 2003) 。 因 此 从 特异 性 位 点 可 以 看 出 ,IR8a 
和 IR25a 更 接近 传统 iGluRs 中 的 non-NMDA 受 体 ， 
建议 将 它们 归 为 non-NMDA 类 型 ,并 命名 为 Putative 
受 体 。 
2.5 iGluRs 和 IRs 三 级 结构 比较 

以 中 华 按 蚊 AsGluRIIa，AsIR2$a，AsIR8a 和 
AsIR60a 为 代表 比较 传统 iGluRs, IR25a, IR8a 和 
IRs 氨基 酸 序 列 的 三 级 结构 (图 4)。AsGluRlla 和 
AsIR25a 氨基 酸 序列 的 三 级 结构 从 N 末端 到 C. 末端 
依次 为 氨基 酸 末 端 区 域 (ATD) 、 配 体 结合 区 (LBD) 
和 离子 通道 区 (IC) 。 其 中 ATD 以 两 边 对 称 的 形式 
组 织 , 两 个 “ 近 端 的 "ATD 以 轴 对 称 线形 成 一 个 小 的 
分 界面 ,两 个 “ 远 端的 ATD "彼此 并 不 接触 ;LBD 是 
由 两 个 部 分 (S1 , S2) 构成 的 “ 捕 晶 草 ” 结构; 而 IC 则 
主要 由 a 螺旋 组 成 , 呈 4 倍 折 三 对 称 。AsIR8a 的 三 
级 结构 除了 有 完整 的 LBD A IC 外 ,其 N 末 端 也 有 
ATD 的 结构 ,只 是 更 分 散 ;而 IR60a 则 只 有 离子 通道 
区 完整 的 IC 和 分 散 的 LBD 区 ,没有 ATD 区 域 。 结 
果 表 明 ,IR8a 与 IR25a 的 三 级 结构 与 传统 的 iGluRs 
更 接近 。 
2.6 iGluRs 和 IRs 的 系统 发 育 关 系 

将 中 华 按 米 冈比亚 按 蚊 、 埃 及 伊 蚊 、 致 倦 库 蚊 
和 黑 腹 果 蝇 5 种 昆虫 的 IR8ga 和 IR25a( 共 10 个 序 
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图 3 5 种 昆虫 代表 性 non-NMDA, NMDA 和 IRs 氨基 酸 序列 的 结构 域 和 特异 性 位 点 注释 
Fig. 3 Domains and specific site annotation of amino acid sequences of representative non-NMDA, NMDA and IRs in five insect species 
A: 5 种 昆虫 代表 性 non-NMDA, NMDA 和 IRs 氨基 酸 序列 多 重 比 对 Multiple alignment of amino acid sequences of representative non-NMDA, NMDA 
and IRs. As; 中 华 按 蚊 Anopheles sinensis; Ag: X] EGET Anopheles gambiae; Dm: Drosophila melanogaster. non-NMDARs, NMDARSs 和 IRs Roue 
称 分 别 用 红色 、 蓝 色 和 黑色 标示 ， 配 体 结合 位 点 有 黑色 框 标示 ; 蓝 色 线条 为 Lig_Chan 结构 域 ,由 跨 膜 区 MI, M2, M3 和 离子 通道 孔 P 组 成 。 黑 
色 、 灰色 和 和 白色 阴影 的 氮 基 酸 序列 其 一 致 性 分 为 别 为 100% , 70% ~99% 和 70% 以 下 。The names of amino add sequences of non-NMDARs, 
NMDARs and IRs are in red, blue and black, respectively, and the black-lined boxes are ligand-gated sites. Lig chan domain is marked with blue lines 


above, consisting of transmembrane region M1 , M2 and M3 and ion channel pore P. The amino acid sequences with 10096 , 70% — 9996 and below 7096 
identity are denoted with black, grey and white shades, respectively. B: 5 种 昆虫 代表 性 non-NMDA, NMDA 和 IRs 结构 域 和 特异 性 位 点 The 
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features of m and specific sites of representative non-NMDA, NMDA and IRs of five insect species. ATD; 44 dE We X tm [X JY, Amino-terminal 
domain; S1, S2; 配 体 结合 区 Ligand-gated sites; Lig Chan; 离子 通道 区 Ion channel region. NMDA 受 体 的 特异 性 位 点 : : PBP2_iGluRs_NMDA, 包 
括 NMDA_ (cd13719/cd06378/cd13720) 3 种 ;non-NMDA 受 体 的 特异 性 位 点 : Kainate 受 体 的 特异 性 位 点 是 PBP2_iGluRs_Kainate 
(cd13714) ,AMPA 受 体 的 特异 性 位 点 是 PBP2_iGluRs_AMPA (cd13715 ) Putative 受 体 的 特异 性 位 点 是 PBP2_iGluRs_Putative (cd13717 ) The 
specific sites for NMDA iGluRs; PBP2 iGluRs NMDA with NMDA_Nr1/Nr2/Nr3 (cd13719/cd06378/cd13720 ) ( three types) ; the specific sites for 
non-NMDA iGluRs: PBP2 . iGluRs, Kainate ( cd13714 ) for Kainate, PBP2. iGluRs, AMPA (cd13715) for AMPA, and PBP2 iGluRs. Putative 


(cd13717) for Putative receptors. 
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图 4 传统 iGluRs, IR25a, IR8a 和 IRs 代表 性 氨基 酸 序列 的 三 维 结构 图 


Fig. 4 Predicted 3D structure of representative amino acid sequences of traditional iGluRs, IR25a, IR8a and IRs 


ATD: 氨基 酸 末 端 区 域 Amino-terminal domain; LBD; 配 体 


列 ) ,代表 性 的 传统 iGluRs (24 个 ) 及 代表 性 IRs 
(30 个 ) 的 氨基 酸 序列 ,分 别 用 最 大 似 然 法 ( ML) , 贝 
叶 斯 (BI) 和 最 大 简约 法 ( MP) 构建 系统 发 育 树 , 结 
果 表 明 这 几 种 方法 得 到 的 结果 基本 一 致 。ML 系统 
发 育 树 分 析 表 明 ( 图 5) ,各 节点 的 自 展 值 (bootstrap 
value) 大 多 在 5096 以 上 ,可 信和 度 高 (Wiens and 
Servedio, 1998) ,与 BI 和 MP 进化 树 一 致 , 黑 腹 果 蝇 
non-NMDA 受 体 中 出 现 了 旁 系 同 源 受 体 ,包括 
AMPA 受 体 中 的 DmClu-R1B 和 DmGluRI, Kainate 
受 体 中 的 DmCG9935 和 DmCG11155, 以 及 
DmGluRIId 和 DmCG5621。 其 他 受 体 在 不 同 物种 中 
则 都 为 直系 同 源 受 体 。IR8a 和 IR25a 与 传统 
iGluRs 的 non-NMDA 受 体 聚 在 一 支 ,它们 与 NMDA 
受 体 构成 一 个 大 的 分 支 。3 种 进化 树 的 不 同 点 在 
于 ,ML 树 与 BI 树 中 ,IRs 形成 了 一 个 单独 的 文系 ， 
而 MP 树 中 的 IRs 基因 则 形成 了 3 SUR. m. 














结合 区 Ligand-binding domain; IC; 离子 通道 区 Ion channel. 


从 ML 树 中 可 以 看 到 ,IR8a 和 IR25a 与 non-NMDA 
fil NMDA 受 体 间 的 遗传 距离 明显 短 于 与 IRs 间 的 距 
离 。 这 些 从 系统 发 育 关系 和 遗传 距离 说 明 , IRSa 和 
IR25a 与 传统 iGluRs 中 的 non-NMDA 受 体 最 接近 。 
2.7 基因 保守 性 和 替换 率 的 比较 分 析 

本 研究 以 中 华 按 蚊 传统 iGluRs AED], IRSa 和 
IR25a 基因 ,以 及 IRs 基因 作为 对 照 , 计 算 了 冈比亚 
按 蚊 、 埃 及 伊 蚊 和 致 伴 库 蚊 中 这 3 类 基因 的 dy/ds 
比值 wa。 如 图 6(B) 所 示 , 传 统 iGluRs 基因 ,IR8a 和 
IR25a 基因 ,以 及 IRs 基因 这 3 类 基因 的 o 值 都 远 
小 于 1, 其 中 传统 iGluRs 基因 与 IR8a 和 IR25a 基因 
的 o 平均 值 接近 ,分 别 为 0.0623 和 0.1038, 说 明 它 
们 在 强 的 纯化 选择 下 进化 ,可 能 编码 功能 性 的 受 体 。 
而 IRs 基因 o 的 平均 值 为 0.2110 ,在 较 弱 的 纯化 选 
择 下 进化 ,因此 可 能 包含 更 多 的 正 向 选择 位 点 ,执行 
不 同 的 功能 。 
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analysis to show the relationships of three types of genes. 


3 讨论 


受 体 (IRs) 是 一 类 高 度 分 化 的 离子 谷 氨 酸 
dun 嗅觉 相关 ,是 化 学 感应 受 体 中 一 个 
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程 中 有 广泛 的 功能 ( Benton, 2007 ) 。 
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中 ,由 于 IR8a 与 IR25a 基因 在 触角 中 表达 ,可 能 与 
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Fig. 6 Comparisons of amino acid length variation ( A) and dy⁄/ds(œ) values (B) in traditional iGluRs, IR8a & IR25a and IRs 


面 与 传统 iGluRs 更 为 相近 ,因此 IR8a. 与 IR25a 的 分 
类 地 位 一 直 是 一 个 未 解决 的 问题 。 本 研究 用 生物 信 
息 学 方法 鉴定 了 4 种 蚊虫 IRSa 与 IR25a 基因 以 及 代 
表 性 的 传统 iGluRs 和 IRs 基因 。 通 过 比较 这 3 类 基 
因 在 氨基 酸 序列 长 度 、 结 构 域 . 配 体 结合 位 点 保守 性 
和 选择 压力 这 4 个 方面 的 特征 可 以 看 出 ,IR8a 和 
IR25a 更 接近 传统 的 iGluRs ; 从 特异 性 位 点 和 系统 发 
育 关系 这 两 个 方面 的 特征 可 以 看 出 ,IR8a 和 IR25a 更 
接近 传统 iGluRs 中 的 non-NMDA 受 体 。 因 此 本 研究 
确定 了 IR8a 和 IR25a 的 特殊 性 ,将 其 重新 定位 为 
non-NMDA 类 型 的 受 体 , 并 命名 为 Putative 受 体 。 

从 氨基 酸 序列 长 度 的 平均 值 看 ,1R8a 和 IR25a 
与 传统 的 iGluRs 明显 更 接近 ,因此 可 能 有 更 相似 的 
功能 。 之 前 的 研究 表明 , iGluRs 的 non-NMDA 
(Kainate 和 AMPA) 和 NMDA 亚 类 在 序列 长 短 上 的 
分 化 会 使 其 不 断 适 应 突 触 交流 机 制 ,从 而 适应 复杂 
的 神经 元 系统 ( Croset et al., 2010) 。 其 次 ,从 结构 
域 和 3D 结构 上 看 ,IR25a 有 与 传统 iGluRs 一 样 负责 
道 的 组 装 以 及 结合 辅助 因子 的 氨基 酸 末 端 区 域 
(ATD) ,IR8a 同样 也 有 类 似 的 区 域 ,只 是 分 散 性 较 
强 , 没 有 预测 到 ;IRs 则 没有 ATD 区 ,这 与 之 前 的 研 
究 相 符合 (Croset et al., 2010)。 所 有 传统 的 iGluRs 
都 具有 3 个 直接 结合 谷 氨 酸 的 氨基 酸 位 点 (R, T 和 
D/E) ,IR8a 同时 具有 这 3 个 位 点 ,IR25a AFAR 
CR) 和 天 冬 氮 酸 (D) ,而 IRs 配 体 结合 位 点 极度 分 
散 ,3 个 位 点 都 没有 。 此 结果 与 之 前 的 研究 结果 一 
致 , 即 黑 腹 果 蝇 、 吏 豆 蚜 、 赤 拟 谷 盗 ` 灰 翅 夜 娥 和 家 和 奉 
IRs 配 体 结合 区 非常 分 散 , 有 可 能 结合 不 同 配 体 , 执 
行 不 同 功 能 ( Benton et al., 2009; 郭 梅 , 2014) 。 同 
时 ,我 们 通过 计算 非 同 义 突变 率 与 同 义 突 变 率 的 比 
值 (dv/ds, w) 来 判断 蛋白 质 编码 基因 的 保守 程度 
以 及 其 是 否 存 在 选择 压力 。 一 般 认为 , 同 义 突变 不 
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受 自 然 选 择 , 而 非 同 义 突变 受到 自然 选择 的 作用 , 若 
w«1 ,表示 所 人 研究 的 基因 属于 纯化 选择 ,wy 之 1 表示 
基因 正 向 选择 , 当 ow —1 时 ,表示 基因 在 进行 中 性 选 
择 ( Swanson and Vacquier, 2002; Swanson et al., 
2003)。 本 研究 结果 发 现 , IR8a 和 IR25a 与 传统 
iGluRs 基因 o 的 平均 值 更 接近 , 远 小 于 1, 它们 在 强 
烈 的 纯化 选择 下 进化 ;而 IRs 基因 o 的 平均 值 较 大 ， 
在 较 弱 的 纯化 选择 下 进化 。 证 明 IRSa 和 IR25a 基 
因 与 传统 iGluRs 基因 一 样 ,保守 性 强 , 可 能 编码 功 
能 性 的 受 体 。 另 外 ,IR8a 和 IR25a 与 传统 iGluRs 都 
属于 胞 质 结合 蛋白 2 超 家 族 的 成 员 , 它 们 的 氨基 酸 
序列 与 1GluRs 有 相似 的 功能 性 特异 位 点 ,其 中 单独 
与 non-NMDA 受 体 相似 的 位 点 有 12 个 ,但 没有 与 
NMDA 受 体 单独 相似 的 位 点 ,因此 IR8a 与 IR25a 更 
接近 传统 iGluRs 中 的 non-NMDA 受 体 。 本 研究 首 
次 利用 了 离子 谷 氨 酸 受 体 氨基 酸 特异 性 位 点 这 个 特 
征 来 确定 基因 之 间 的 相似 性 ,并 根据 位 点 相似 性 来 
进行 分 类 ,这 为 了 8a 与 IR25a 的 重新 定位 提供 了 直 
接 证 据 。 系 统 发 育 分 析 也 表明 ,IR8a 和 IR25a 与 
non-NMDA 聚 在 一 文 上 ,遗传 距离 也 与 其 相近 , 远 小 
于 IRs。 这 与 之 前 黑 腹 果 蝇 (Bengtsson et al., 2012; 
Liu et al., 2012) ) 版 虫 和 棉铃 虫 (Liu et al., 2012; 
Glaser et al., 2013) 等 物种 离子 谷 氨 酸 受 体 的 系统 
发 育 关系 结果 一 致 。 

本 人 研究 从 受 体 的 结构 域 、 特 异性 位 点 和 系统 发 
育 关系 等 各 个 方面 比较 结果 表明 ,IR8a 和 IR25a 与 
传统 iGluRs 更 为 接近 ,因此 可 能 有 与 其 一 样 在 神经 
元 中 发 挥 作用 的 功能 ,这 可 为 以 后 的 相关 实验 验证 
工作 提供 假设 。 本 研究 丰富 了 蚊虫 离子 谷 氮 酸 受 体 
的 基因 数据 ,并 有 利于 推动 离子 谷 氮 酸 受 体 基 因 在 
蚊虫 生长 发 育 中 基因 的 表达 调控、 功能 及 其 他 生物 
学 方面 的 研究 。 
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